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Avertissement : Nous introduisons 
l’électromagnétisme dans ce cours en nous basant 
sur les deux enseignements précédents de L1 et de L2 au 
cours desquels ont été développés l’électrostatique et la 
magnétostatique, théories que nous supposerons acquises 
(c'est-à dire que les équations V -E = p/€0, V x B = uoj, 
VxE=0et VB = 0 doivent être connues et leurs 
conséquences comprises). Nous verrons ainsi les modifica- 
tions qu’il faut apporter aux équations de Maxwell pour LAN 
prendre en compte de nouveaux phénomènes essentiels 
liés aux aspects dépendant du temps. 





Sous la statue de Maxwell, 
à Edimbourg. 


1. Introduction 


Cette section introductive a pour but de donner aux lectrices et aux lecteurs une idée de la 
richesse de l’électromagnétisme sur le plan mathématique par la diversité des formulations que 
l’on peut en donner. Il ne s’agit pas à ce niveau de comprendre ! Certaines des notations utilisées 
dans cette section vous sont totalement étrangères. Mais comme pour l’apprentissage d’une langue 
étrangère, la physique demande parfois une plongée en immersion pour donner l'envie d’aller plus 
loin... 

e L’électromagnétisme est une théoie physique exemplaire à bien des égards. C’est la première 
fois que sont unifés des phénomènes antérieurement supposés indépendants les uns des autres (les 
domaines de Pélectricité et du magnétisme). Elle sert à ce titre de premier exemple d’unification 
entre les domaines liés à l’existence de E et ceux liés à l’existence de B : 


E 
VE: V x B = poj + Loco (1) 
VxE--Æ, V-B=0. (2) 
avec en conséquence 
CA 
E = -Vo - — 3 
Vo, (3) 
B=Vx%xA, (4) 
Op 
ss +1 = (0. 5 
Sd | (5) 
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Tous les champs qui apparaissent dans ces équations dépendent dans le cas général des coordonnées 
spatiales et du temps. Nous verrons tout cela ensemble. 

e C’est une théorie relativiste (au sens de l’invariance par changements de référentiels galiléens 
de la relativité restreinte). Les équations de Maxwell s’écrivent alors 


OuE = uoj”, Fur = OuAr — 0yAu, An = MuvA” (6) 
A = (b/c, —AÀ), Ju . (pc, —i) (7) 


avec pour conséquence 0,3" = 0. Vous verrez cela en L3. Sous forme développée cela conduirait 
plutôt à écrire les équations (1) et (2) comme 
1 0E 
: =; 12> ir = 1 
V:E=r, VxB- = Hoi (8) 
VXxE+Æ-=-0, V-B=0, (9) 


e C’est une théorie de jauge et de ce point de vue elle sert de précurseur aux théories de 
jauges des interactions fondamentales dans le modèle standard. Cela signifie que dès lors que 
les propriétés de la matière peuvent être décrites au moyen de champs complexes (mécanique 
quantique), l’invariance par changement de phase est une symétrie globale de la théorie qui induit 
une loi de conservation, celle de la charge électrique et que l’extension de cette symétrie aux 
changements de phase locaux permet de construire une action covariante de jauge et invariante 
relativiste qui engendre l’électromagnétisme de Maxwell : 


1 
— Jéz (- nr" - At) = 0: (10) 





Vous verrez cela en M1. 
e C’est un théorie géométrique et topologique qui se prête à des formalismes mathématiques 
modernes, comme celui des formes différentielles : 


œp+dj=0, (11) 
2 3 

D = bp; (12) 
1 2 2 

dH = J+@D, (13) 

2 1 
&B+dE = 0, (14) 

2 1 
D = Ee(x3E), (15) 
1 2 

H —= ul(x3B). (16) 


Vous verrez peut-être cela un jour en détail, mais comme cadeau de fin d’année je vous offrirai une 
introduction à ce formalisme très élégant en fin de cours ! 


2. Notations 


1.- Les points de l’espace R° sont génériquement notés comme des vecteurs à 3 dimensions, 
en caractères gras, x. La variable de temps est notée t. 

2.- Les champs scalaires (potentiel scalaire) sont notés d(x, t) et les champs vectoriels A(x,t), 
B(x,t) ou E(x, t). Les sources des champs sont également des champs, scalaire p(x, t) ou vectoriel 
ic. 

3.- Les opérateurs vectoriels usuels sont notés comme grad 4 = Vo, Ag = V°6, div A = VA, 
rot A = V x A et AA = V°A. On utilise fréquemment les propriétés utiles V?9 = div(grad d) 
et rot(rot À) = grad(div A) — V’A. 
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4. - Théorème de Green-Ostrogradski, intégrale d’une divergence sur un volume, 


| divE av fB.as. (17) 
@ 
on 


Q est un volume de R* et 0Q son bord (une surface fermée) orienté vers l'extérieur. Le second 
membre ci-dessus est un flux, noté ®;,Q(E) (la notation ® (majuscule) ne doit pas être confondue 
avec celle du potentiel scalaire, @ !). 

5.- Théorème de Green-Stokes, flux d’un rotationnel à travers une surface, 


Jrots.us- {8.4 (18) 
* 0x 
X est une surface orientée plongée dans R* et 02 son bord (un contour fermé orienté également, 


la convention d'orientation de la surface et du bord est la convention de compatibilité standard). 
6.- Circulation d’un gradient 





| grad 4 : dé = Ÿ ao = (Az) — (Ai). (19) 
Â 


0A 


À est une ligne orientée plongé dans R* et OA son bord (deux points ordonnés). 


3. Les théorèmes généraux de l’analyse vectorielle et les opérateurs différentiels 
associés 


3.1. Notion de champ 


Soit R° l’espace euclidien ordinaire que l’on appelera M. Les points de M sont repérés par 
x € M. Une fonction à valeurs réelles définie en tout point de M, f(x) ER, est un champ scalaire 
réel. Un vecteur A(x) défini en chaque point de M par ses composantes dans une base (qui sont 
autant de champs scalaires réels) définit un champ vectoriel 


A(x) = | 4,(x) (20) 


avec ici le choix d’une base cartésienne unitaire {u,,u,,,u,} pour décomposer A(x). 
Exercice 1 : Coordonnées cylindriques et sphériques. 


Ecrivez également la décomposition en coordonnées cylindriques et sphériques. Faites soigneuse- 
ment le schéma de ces systèmes de coordonnées. 


3.2. Comment un champ scalaire varie-t-il ou la notion de gradient 


Soit f(xo) un champ scalaire réel et x0 € M. En xo, le champ scalaire prend la valeur f(xo) et 
on cherche à savoir ce qu’il vaut en un point x voisin de xo. Généralisant la notion de dérivée d’une 
fonction d’une variable, c’est la notion de gradient qui permet d’écrire le développement limité au 
premier ordre : 


f(x) = f(&o) + (x —x0): grad f|,, + O(Ix — xo|”). (21) 


grad f, ou Vf, est un champ de vecteurs, ici évalué au point xo. Si x et x9 sont infiniment proches, 
x = xp + d£, et si x est sur la même surface iso-valeur de f que xo, alors on en déduit 


de - grad f|,, —0, (22) 








ce qui montre que grad f Le est perpendiculaire à cette surface iso-valeur en x9. Sinon, dans le cas 
général, f varie entre x et xo de 
JC) = f(&o) + df (23) 
avec df la différentielle de f en x9 de sorte 
df(&o) = dé: grad f |, . (24) 


Cette relation permet de calculer facilement grad f dans les systèmes de coordonnées orthonor- 
maux usuels. Soit (q1,Qq2,q3) un tel système de coordonnées orthogonales rapportées à une base de 
vecteurs normés et orthogonaux entre eux, (u:,u,us), u; -u; = 6;; 


X = Qi + QU + q3U3, dl = dlyu:1 + dlous + d{;u3 (25) 
(les dl; coïncident avec les dg; en cartésiennes mais en diffèrent en général). On a donc 


3 


F7 ja 
i=1 
3 
nn. i grad f|,,) (26) 


soit en tout point 


3 

1 of 
grad f = — —— 
2 dl; dqi 


Par exemple en coordonnées sphériques (r, 0,4), 





dl = dru, + rdôus + r sin dou, (28) 
ce qui conduit à 
1of,. 1 Of 1 of 
df — 2 Qu + ——— "4 
grad f on T0 F0 do 
Of 10f 1 of 
= # s 2 
&r Tao MAN TE D @) 


On vérifie aisément la dimensionnalité de chacune des composantes, homogènes à [dim f|/[longueur] 


Exercice 2 : Gradient en coordonnées cylindriques. 
Faites le calcul en cylindriques. 


3.3. Comment un champ de vecteurs traverse-t-il une surface ou la notion de divergence 


On définit tout d’abord le flux d’un champ de vecteurs A(x) à travers une surface Y. Locale- 
ment, sur une portion infinitésimale orientée dS de Y, le champ de vecteurs peut être considéré 
comme constant et on peut évaluer le produit scalaire 


A(x) - 48 = |A(x)| || d8 || cos(A(x), 4). (30) 


Cette quantité est nulle si le champ de vecteurs est localement parallèle à la surface et sera d’autant 
plus grande (en valeur absolue) que A(x) aura une composante importante perpendiculairement 
à cette surface. Le flux total du champ de vecteurs A à travers la surface Z, Ss(A), s'obtient par 
intégration, 


A) = Î A(x) : dS. (31) 


Il mesure la manière dont le champ de vecteurs traverse la surface. La surface est une grandeur 
orientée, le flux sera positif ou négatif selon la convention choisie pour orienter X. Dans le cas où 
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la surface est le bord d’un volume Y = 00, on parle de surface fermée et la convention usuelle 
est d'orienter les dS vers l’extérieur du volume (. On définit dans ce cas la divergence de A, 
div A(x) ou V - A(x) pour mesurer avec quelle intensité A traverse localement une surface fermée 
0Q limitant un volume dV entourant x, 


1 
div A(x) = En Pin } A(x) - dS (32) 
eN(x) 


De manière explicite en coordonnées, on a 


1 
div A(x) = lim A 
TE mer GE 1dtzdts)dq + 





© © 
n Sat PAIE e -(Aadhdtz)dm ). 


(33) 
par exemple en coordonnées sphériques on aurait qg1 = r, g2 = 0, gs = vw, dl = dr, dl = rdb, 
dl3 = r sin 0dç et 








1 
iv À =. fn 
ANA di 0 drrdôr sin dy 


(£ (A,rdôr sin 0dçp)dr + 2 (4er sin Üdçpdr)d0 + = (4e drrdbhds). (34) 


soit finalement 


ol 1 à 1 à 
Eau Ar) En r sin 0 00 PIE LL DS 





div A(x) = : (A). (35) 


r sin Ô êg 
La divergence d’un champ de vecteurs est d’autant plus intense que le champ de vecteurs 
traverse les surface d’un volume fermé. 


Exercice 3 : Divergence en coordonnées cartésiennes et cylindriques. 
À titre d'exercice il est recommandé de traiter les cas des coordonnées catésiennes et cylin- 
driques. 


3.4. Comment un champ de vecteurs circule-t-il le long d’un contour fermé ou la notion de 
rotationnel 


La circulation d’un champ de vecteurs le long d’un contour fermé infinitésimal orienté qui limite 
une surface Z et supposé être constitué de branches alignées le long des éléments de longueur de 
base (d’une base orthogonale pour ce qui nous intéresse ici) s'écrit 





An A(x)-d£ =  (Ao(q + dqi,q2,q3)dl2(qi + dq1, 2, q3) — A2(q1, q2, q3)dl2(q1, q2, q3)) 
D 
—(A2(q1,q2 + dq2, q3)dl1(q1, 42 + dq2, q3) — A2(q1, q2, 43)dl1(q1, q2, q3)) 
© 0) 
= {hd (Aide), (36) 
og1 Oq2 


Attention, il y a une convention d'orientation qui est indispensable et sous-entendue dans l’expression 
ci-dessus qui pourrait prendre un signe opposé. 

Par définition, cette circulation est égale au flux du champ de vecteurs rot A(x), ou V x A(x), et 
dans l’exemple ci-dessus où la base est orthogonale, ce flux est simplement donné par la composante 
3 du rotationnel car la normale à Y est portée par l’axe 3, soit 


(rot A(x) : us )d£: dla = im A(x) $ dE, 
| 02 
1 Ô Ô 
t A(x) : = —— A A ; 
GotAG)w) = gr (ar (Ardie) - (Ardh)) 7) 














Appliquée au cas des coordonnées sphériques avec g1 = r, g2 = 0, g3 = w, dl: = dr, dl = rdb, 
dl3 = r sin 0d, on obtient 

10 10 
EC VS PE 
r Ôr (en r 00 


On comprend que le rotationnel d’un champ de vecteurs est d'autant plus intense que le champ 
de vecteurs tourne quand on se déplace le long du contour fermé. 


(rot A(x) - u,) = Ar. (38) 


Exercice 4 : Les autres composantes de rotationnel. 
Déterminez les composantes de rotationnel suivant u, et ugs en coordonnées sphériques, puis 
les trois composantes en cartésiennes et en cylindriques. 


4. L'induction électromagnétique 


4.1. les faits expérimentaux 


La situation caractéristique mettant en évidence le phénomène d’induction électromagnétique 
est celle où un permier circuit (circuit 1) actif est parcouru par un courant électrique et crée un 
champ magnétique dans l’espace. Un second circuit (circuit 2), passif, c’est-à-dire dépourvu de 
générateur se voit néanmoins, dans certaines situations que l’on va décrire, également siège d’un 
courant électrique, mesuré par exemple au moyen d’un galvanomètre, qu’on appelle courant induit. 

Dans une expérience 1, on déplace ou on déforme le circuit 2 dans le champ du circuit 1. Un 
courant induit apparaît dans le circuit 2. 

Dans l’expérience 2, le circuit 2 est immobile, mais c’est le circuit 1 qui est déplacé. Cette 
seconde expérience peut correspondre à la même que la première, mais vue d’un autre référentiel 
et de ce point de vue elle n’apporte pas de situation physique réellement différente. 

Dans l’expérienc 3, les circuits 1 et 2 restent immobiles, mais le circuit 1 est parcouru par un 
courant 1(t) variable et crée donc un champ magnétique B(t) variable dans le temps. 

Dans l’expérience 1, on a @,B = 0 alors que dans les expériences 2 et 3, 4 B + 0. Dans tous 
les cas il apparît nécessairement une f.e.m d’induction que nous allons définir pour qu’un courant 
circule dans le circuit 2. 





4.2. Notion de force électromotrice d’induction 


Considérons un circuit fermé l'. Un électron dans ce circuit ressent une force motrice Fu, 
puisqu'il est mis en mouvement. En effet, dans un circuit en régime permanent, les charges sont 
en mouvement à vitesse constante à condition qu’une force motrice agisse sur eux, de sorte que 
l'équation de mouvement incluant l’effet de la résistance du conducteur, 





du(t) m 
= Fu — t 39 
ME = Fu — u(b, (39) 
admette comme solution stationnaire u = T Fu. Par définition, la f.e.m (l’appellation de force est 
un peu abusive, puisqu'il s’agit en fait d’une différence de potentiel (d.d.p), mais c’est l’usage,) est 
la d.d.p correspondant à la force motrice qui déplace les électrons, 


F F 
e = } Se de, A] = champ électrique. (40) 
q q 


e mesure ainsi l’énergie reçue par le biais de Fu par une charge unité faisant le tour du circuit. 
Dans un circuit actif, c’est le générateur qui fournit la force motrice. Elle peut être d’origine 
chimique (batterie), lumineuse (capteur solaire), mécanique (piezo-électrique), … Dans le circuit 
passif, la seule cause possible à la force motrice et donc à la f.e.m d’induction est le champ élec- 
tromagnétique qui règne au voisinage du circuit 2, ce ne peut donc être que la force de Lorentz, 
Fu 


=E+vxB. (41) 
q 
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Pour traiter les différents cas de figure, la vitesse des porteurs de charge peut être décomposée en 
V=vr+u (42) 


où vr est la vitesse de l’élément de circuit 2 par rapport au circuit 1 et u est la vitesse des électrons 
dans le circuit. Comme u/||d£, il reste uniquement 


Pr $e PR (43) 
r 


4.3. Les lois de Faraday de l'induction 
4.3.1. Circuit passif mobile dans un champ B constant 


C’est la cas de l’expérience 1 avec 0,B = 0. Si l’on suppose, comme en électrostatique, que 
rot E = 0, alors 


je -d=0, et ejna = 10 X B) - d£. (44) 
r r 


On verra plus tard, dans le formalisme plus général, que l’hypothèse sur le rotationnel de E est 
correcte car nous sommes en régime ©B = 0. Or vr peut être écrite comme la variation de la 
position d’un point M du circuit au cours de son déplacement infinitésimal, vr = JON, de sorte 


que 


nee OM ne OM gp OMS (5) 


et dOM x d£ = ndS où dS est la surface balayée par d£ au cours du déplacement infinitésimal 
du circuit et n est la normale orientée à cette surface (dS est bien un infiniment petit du second 
ordre, mais on n'écrit pas d?$S en général, comme on n’écrit pas d’V pour un élément de volume). 
Il vient donc une intégrale de surface, étendue à la surface Y limitée par le contour F, L = ©, 


_dbe 
dt 





Gna = —ÿ; | B-nds = (46) 
où l’on à défini ®. le flux du champ magnétique coupé par le circuit au cours de son déplacement. 
La variation he n’est rien d’autre que la variation du flux a à travers le circuit entre l’état avant 
déplacement ou déformmation et l’état après. 

On note qu’il n’y a de définition non ambiguë de Y qu'avec le flux coupé (Y est alors la surface 
balayée par le circuit au cours de son déplacement), mais en raison de la propriété divB = O, 
l'intégrale pour calculer ® est indépendante du choix de la surface qui s’appuie sur let la définition 
de e est sans équivoque. 





4.3.2. Circuit fixe dans un champ B variable 


Cette fois vr = 0 et la f.e.m doit se réduire à 


Gina = PE de +0 (47) 
T 





ce qui entraîne que rot E + 0 et donc que E ne dérive plus d’un potentiel scalaire. Si l’on met en 
regard le fait que dans cette expérience rot E + 0 et ;B + 0 alors que l’on avait dans la précédente 
rot E = 0 et ©,B = 0 et si l’on suppose une relation entre ces deux faits on peut poser de manière 
générale une relation de la forme 











rotE fonction de 0,B 
= _f(0) + ac;B + termes d’ordre 2 


= a0;B (48) 








car on à f(0) = 0 pour satisfaire à la situation de l’expérience 1 et si l’on retient pour l’électro- 
magnétisme une théorie linéaire comme l’expérience générale le confirme, les termes d’ordre 2 et 
au-delà doivent s’annuler également. 

Supposons cette relation satisfaite, que pouvons-nous en déduire ? 


Eind — few | rot E : dS 
: E(T) 


a | B : dS 
E(T) 


d® 


ce qui demande que a = —1 pour rendre compte des résultats de l’expérience 1. 
On peut donc écrire les lois de Faraday sous forme locale ou intégrale : 


d® 


rotE = —-—, où eng = ——. 
ind dt 


É 0) 


On en déduit que la relation entre les champs et les potentiels prend une forme plus générale 
également. La relation div B = 0 demande que 


B = rot A (51) 
et rot E = -2 s'écrit alors rot(E + 0, A) = 0 ce qui demande que E + 0,A = —V, soit 
E = -Vo — +0,A. (52) 


Cette nouvelle forme du champ électrique en fonction des potentiels décrit à la fois les phénomènes 
de l’électrostatique et ceux qui dépendent du temps. 


4.4. Courants de Foucault 


Nous allons traiter une application aux courants de Foucault dans le cadre d’une approximation 
que nous améliorerons une fois démontrée la forme complète de la dernière équation de Maxwell 
dans la section suivante. On considère donc un cylindre de rayon a de longueur infinie (ou bien de 
longueur finie mais on ne tient pas compte des effets de bords, de sorte que la symétrie cylindrique 
est satisfaite avec invariance par translation le long de l’axe Oz du cylindre et par rotation d’angle 
0 autour de cet axe). Ce cylindre est soumis à un champ magnétique externe variable qui s’écrit 
en première approximation 

B(t) = Bocoswtu, = B,(t)u.. (53) 


On verra plus tard que cette forme n’est pas compatible avec l’ensemble des équations de Maxwell, 
mais la symétrie demande au moins une forme telle que 


B(r,t) = B,(r,t)u. (54) 


On s'intéresse au champ électrique induit par ce champ magnétique et pour cela on utilise la loi 
de Faraday 
VXxE = -0,B (55) 


avec E(x) = Eo(r,t)us en raison des symétries du problème. Au lieu de la forme complète de 
l’opérateur rotationnel en coordonnées cylindriques, il est avantageux ici, en raison des symétries 
des champs, d'utiliser le théorème de Stokes, 


Jvxe.as-fe.æ--f anus. (56) 
S S 


os 
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Pour exprimer facilement le dernier terme, on choisit pour surface l’élément dS = dr rd@ u, 
dont on oriente le contour en conséquence. Il vient pour le flux élémentaire de —0,B l’expression 
— 04 Bo coswt rdrdû = Bow sinwt rdrd@. La circulation de E le long du bord de cette surface élé- 
mentaire vaut (r + dr)d0 Eo(r + dr,t) — rd0 Eo(r,t) = À (rEo(r,t)) drd0. L'égalité conduit à 


2 (rEo(r,t)) = Bowr snwt (57) 
ñ 


que l’on intègre en 


1 
Eo(r,t) = 3 Bowr sinwt. (58) 


Ce champ électrique à l’intérieur du conducteur de conductivité o donne naissance à des courants 
électriques appelés courants induits ou courants de Foucault, 


1 
j = jo(r,t)uo = D Pour sin wtUy (59) 


qui dissipent une puissance j : E = ÀBiow?r? sin” wt par effet Joule. En moyenne sur le temps et 
intégrée sur le volume du cylindre, cette puissance vaut 


h a 
1 
| de | 27rrdr j + Eh = The" Bou. (60) 
0 0 


4.5. Exercices 


Exercice 5 : Frein TELMA. 

Un cylindre conducteur plein, amagnétique, tourne lentement à vitesse angulaire w constante 
autour de son axe. Il est placé dans un champ magnétique uniforme et constant Bo perpendiculaire 
à l’axe. 

1.- Calculer la densité de courant induit (courant de Foucault). 

2.- Calculer la puissance dissipée par effet Joule par unité de longueur et déterminer le moment 
du couple F. 


Exercice 6 : Coefficient d’inductance propre d’une bobine. 

On considère une bobine assimilée à un solénoïide infini de n spires par unité de longueur (ou 
à un solénoïde fini de N spires au total mais pour lequel on néglige les effets de bords). 

1.- Calculer le coefficient d’autoinductance de la bobine par les méthodes usuelles en magné- 
tostatique : 

— à partir du calcul du flux propre 6, 

— à partir du calcul de l’énergie magnétique U,, emmagasinée par la bobine. 

2.- Calculer enfin ce coefficient en appliquant les lois de l’induction, #.e. à partir de l’expression 
de la f.e.m. d’autoinduction e4g aux bornes de la bobine, déduite elle-même du champ électro- 
moteur E,,. 


Exercice 7 : Courants de Foucault. 

Un cylindre métallique de conductivité ao, d'épaisseur h et de rayon a est soumis à un champ 
magnétique variable de la forme B = B5 coswt où Bo est parallèle à l’axe Oz du cylindre. 

1.- On suppose que w est suffisamment faible pour que l’on puisse considérer valide l’hypothèse 
des régimes quasi-permanents, et considérer que B est uniforme à chaque instant. Calculer la densité 
de courants j induits dans le cylindre. 

2.- Calculer ensuite la puissance dissipée par effet Joule. Evaluer l’intensité qui circule dans 
le cylindre et en déduire la résistance du cylindre. 

3. - La fréquence étant toujours supposée faible, peut-on évaluer la correction au champ 
magnétique due aux courants de Foucault ? (il faut pour cela connaître la fomre généralisée du 
théorème d'Ampère, établie à la section suivante). 








Exercice 8 : Freinage par induction. 

Un cadre rectangulaire peut glisser verticalement sans frottement dans des rails. Il est constitué 
d’un fil de résistance À. Il tombe verticalement et atteint une région de l’espace où règne un champ 
magnétique uniforme et constant, perpendiculaire au plan du cadre. Déterminer la vitesse du 
cadre avant qu'il ne plonge dans le champ, pendant puis après. On néglige dans un premier temps 
le champ B: dû au courant induit dans le cadre. Peut-on ensuite en tenir compte ? 


Exercice 9 : Circuit rectangulaire en expansion. 

On considère le circuit en expansion constitué d’un conducteur mobile se déplaçant sur des rails 
sans frottement, le tout étant placé dans un champ magnétique uniforme et constant. 

1.- Calculer le courant induit si l’on néglige la résistance des rails. 

2.- La résistance des rails étant donnée en (Q/m, soit r, calculer de nouveau l'intensité i. On 
considère qu’à t = 0 le conducteur mobile était en contact avec AB. 

3. - On considère maintenant la situation suivante où la barre a une masse m et où à t = O, 
lorsque l’on ferme l'interrupteur, la barre est encore immobile. 

— Déterminer les variations de vitesse de la barre. Donner la vitesse limite. 

— Cette vitesse limite est atteinte au bout d’un temps {1 pris comme nouvelle origine des temps. 
La barre est alors soumise à une force de frottement constante f. Déterminer v(t). 

—Même question si la force de frottement peut être exprimée par f = —kv. 


5. Généralisation du théorème d'Ampère 


5.1. Le problème de la conservation de la charge 


Au point où nous en sommes, la situation de l’électromagnétisme est décrite par les équations 

rotE=-%Æ, divB=0 (61) 

. + = . 

div E = Ps rot B = Lo) (62) 
en fonction des sources de charges et de courants p et j. Tous les champs sont a priori des fonctions 
de l’espace et du temps. On sait par ailleurs empiriquement que la charge électrique se conserve, 
donc que ces sources obéissent à l’équation locale 


p + divj = 0 (63) 


car si l’on intègre sur un volume Q 


d 
J'apav- | dvjdv = pi-us (64) 
a dt a 

on 


et que l’on étend le volume Q — æ, &,, 3j: dS — zero (il n’y a pas de sources à l'infini) 
et l’on retrouve bien le fait que la charge totale dans l'Univers est conservée. La forme locale 
(sans étendre le volume d'intégration à l'infini) exprime plutôt la conservation dans les processus 
à toutes les échelles : la charge augmente dans un volume donné s’il y a un flux de courant entrant 
et réciproquement elle diminue si le flux net est sortant. 
Or cette équation locale est incompatible dans le cas général avec la forme actuelle du théorème 
d'Ampère qui conduit plutôt à 
div rot B = 0 = yo divij. (65) 


Il est indispensable pour rendre compte de la conservation de la charge électrique de généraliser 
le théorème d'Ampère dans les cas des phénomènes dépendant du temps, tout en préservant sa 
forme dans la cas statique. Pour cela, on est amené à postuler une expression de la forme 


rot B = o(j + B) (66) 





10 


L2 Licence Physique Chimie 





de sorte que divrotB = 0 entraîne div(j + B) = 0. Pour que cette dernière égalité redonne 
la conservation de la charge, il faut que div B = &@p. Si l’on combine ce résultat à la première 
équation de Maxwell, Gp = @(E0 div E), on obtient que 


B = &oûùE (67) 


ce qui conduit à l’expression du théorème d'Ampère généralisé, 
0E 


3 (68) 


rot B = Loj + Echo 


. . à . À . . 
qui revient à ajouter un terme DE au courant de conduction j. Ce n’est pas un courant de 


charges et la dénomination usuelle est celle de courant de déplacement. Ce terme donne naissance 
aux phénomènes de propagation et il est justifié par ses conséquences. Il est remarquable que 
Maxwell ait eu une telle intuition, confirmée plus tard par la découverte expérimentale des ondes 
électromagnétiques par Hertz. 

Les quatre équations de Maxwell sont en définitive les suivantes : 


°B 


rot E = gr Loi de l'induction de Faraday (69) 
divB=0, Conservation du flux magnétique (70) 
divE = Len Théorème de Gauss (71) 
rot B = Jo ( + 2%) ,; Théorème d'Ampère généralisé (72) 


5.2. L’équation des ondes 


La nouveauté essentielle de l’extention du théorème d'Ampère est qu’il permet d'établir pour 
les champs E et B dans le vide et même en l’absence de sources des équations de propagation. En 
effet, en prenant rot rot E on forme 


rotrotE = —@,(u0c0@4E) 
— — 106007 E, (73) 


soit, en utilisant également la propriété rot rot E = grad(div E) — V'E, 


1®E 


vhs 
c2 ot? 


(74) 
De manière parfaitement équivalente, on montre que la même équation de propagation vaut pour 
le champ magnétique si l’on part de rot rot B, 


2 
vB-108 5 (75) 


| € dt 


En présence de sources de charges et de courants, les équations font intervenir un second membre 
(équations inhomogènes) et l’on obtient 


9 1®@E 1 0 
V TT nd (76) 
1 ©B : 








5.3. Les équations pour les potentiels 
5.3.1. Invariance de jauge 


Revenons tout d’abord sur la notion d’invariance de jauge en électromagnétisme. Le point clé 
est ici qu’une situation physique donnée, définie en termes de champs électrique et magnétique qui 
déterminent les forces agissant sur des charges test, peut être décrite par une infinité de manières 
distinctes par des jeux de potentiels différents. 

La conservation du flux magnétique est automatiquement réalisée on le sait si le champ magné- 
tique s’écrit comme le rotationnel d’un champ de vecteurs appelé le potentiel vecteur, B = rot A. 
Mais cela ne suffit pas à définir À de manière univoque car le choix 


A'=A+VX, rotA =B (78) 


définit le même champ magnétique (si y(x,t) est une fonction définie de manière univoque sur 
l'espace et le temps). Ce caractère arbitraire dans la définition de A a suscité bien des interrogations 
et les physiciens au cours du temps ont adopté des positions variées, allant pour Heaviside par 
exemple, d’un reniement de telles quantités monstrueuses, à des positions pragmatiques, comme 
Maxwell, pour qui les potentiels sont des intermédiaires de calcul, la physique étant gouvernée par 
les champs qui déterminent les forces, E et B et donc la dynamique. On peut alors s'interroger 
sur une expression telle que celle du potentiel vecteur d’une distribution stationnaire de courants 


étudiée en magnétostatique, 
H0 JG) dV 
A(x) = : 
Go) AT Î [x — x] (79) 





Cette expression vaut dans une jauge donnée, c’est celle du potentiel vecteur ayant la symétrie 
de la distribution de charges, mais une infinité d’autres choix sont possibles (au-delà de simples 
constantes à l'infini). 

La discussion ci-dessus s’étend au potentiel scalaire ®, lui aussi défini de façon non univoque. 
En effet, si un choix de potentiels 9 et A définissent le champ électrique, E = —-V@— @,A, on peut 
également réaliser ce même champ électrique avec la combinaison 





A'=A+VXx, et #=p-@x -Vé -à@A'=E. (80) 


Le potentiel scalaire d n’est donc également défini qu’à une fonction près dans le cas général et 
une expression sans ambiguïté n’est possible qu’en fixant certaines conditions. L'expression q@ 
ne représente alors l’énergie d’une charge q dans un champ que si d possède les symétries de la 
distribution de charges et s’annule à l’infini pour une distribution finie. 

On parle de liberté de jauge pour exprimer les diverses possibilités offertes par le choix de la 
fonction x et lorsqu'un choix est fait on parle de choix de jauge. Ces choix sont très intimement 
liés aux symétries exploitées dans le problème physique comme on va le voir. L’invariance de jauge 
est une symétrie que l’on impose aux théories physiques pour que les quantités physiques soient 
indépendantes du choix de jauge. Cette approche, initiée en 1918 par le mathématicien Weyl, puis 
popularisée par les physiciens Yang et Mills dans les années 1950 s’est révélée d’une redoutable 
fécondité pour élaborer le modèle standard des interactions fondamentales. 

Le statut physique des potentiels lui aussi à été interrogé longuement et Aharonov et Bohm et 
ont proposé en 1959 des expériences permettant de mettre en évidence la signification physique, 
non pas des potentiels eux-mêmes, mais de certaines quantités intégrales de ces champs. Cela a 
donné naissance à des quantités topologiques invariantes qui jouent de nos jours un rôle central 
dans de nombreux domaines de la physique. 


5.3.2. Une illustration élémentaire et instructive de la liberté de jauge dans un cas d’école 


Je vous propose, chères lectrices et chers lecteurs, une digression par rapport au fil de ce cours 
pour illustrer cette liberté de jauge et les conséquences en termes de symétrie et de quantités con- 
servées dans un cas élémentaire que vous connaissez très bien par ailleurs. Il s’agit du condensateur 
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plan entre les armatures duquel on lâche depuis l’une des armatures et sans vitesse initiale une 
charge q et l’on demande la vitesse acquise lorsque la charge atteint la seconde armature. 

On considère ainsi un condensateur plan parfait composé d’une armature en x = 0 chargée 
négativement et d’une armature plane parallèle chargée positivement et placée en x = d. Il règne 
entre les armatures un champ électrique E = Eu, avec E = U/d où U est la ddp imposée par un 
générateur. La charge q est supposée positive. Le problème est résolu classiquement par le principe 
fondamental de la dynamique que l’on écrit 


du(t) 


: 1 E 2 
dt 


m = 
2m 





=qE, soit Ho 8 a(t) (81) 





où l’on a tenu compte des conditions initiales x(t = 0) = v(t = 0) = 0. On en déduit le temps pour 


atteindre la seconde armature, t4 = 1/ ns et donc la vitesse à cet instant 


gE |2md 2qEd 
ARE SEAITE = A) (82) 


On préfère souvent la solution plus rapide au moyen de la conservation de l’énergie mécanique 
qui demande que la variation d'énergie cinétique K(d) — K(0) = mv?(d) soit compensée par la 


variation d'énergie potentielle V(0) — V(d) = qU, soit 


v(d) = ee (83) 


Cette solution élegante sous-entend certains présupposés que nous allons mettre en valeur grâce 
au formalisme lagrangien. Rappelons que le lagrangien L(x, v,t) est une fonction qui permet de 
décrire la dynamique d’une particule au moyen des équations d’Euler-Lagrange, 








ÔL dôL 
ST (84) 
0x  dtOov 

Ces équations sont équivalentes à la relation fondamentale de la dynamique. Pour un système 
conservatif le lagrangien s'écrit généralement comme la différence entre l’énergie cinétique K et 
l'énergie potentielle V. Pour une particule dans un champ électromagnétique, elle s’écrit cependant! 


Le, v,t) = Smiv(DP — q(p(x) — A(x) : v()). (85) 


Par transformation de Legendre on construit une autre fonction, l’hamiltonien, H(x, p,t) : 


H(x, pt) = p-v— L(x, vi). (86) 
La contrainte _. — 0 entraîne 
| êL 
P — FE (87) 


et la dynamique, équivalent aux équations d’Euler-Lagrange, est alors donnée par 


CH dp CH dx 





= = ; 88 
ox dt’ op dt (88) 

Pour une particule dans un champ électromagnétique, l’impulsion est donnée par 
p = Mmv + GA. (89) 


1C’est un cas un peu particulier car aucune énergie potentielle n’est associée à la force magnétique qui ne 
tavaille pas. On peut se convaincre toutefois que les équations d’Euler-Lagrange redonnent bien avec ce lagrangien 
la dynamique d’une charge soumise à la force de Lorentz. 
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Du point de vue de l’étude des lois de conservation, le formalisme lagrangien est très utile. 
En effet, si l’une des variables de coordonnée q; n'apparaît pas explicitement dans le lagrangien 
L(di, di, t); 9 d 

L D; 
—— =0 alors —Z2=0, ou p; = const 90 
°q dt Pj (90) 
êL 


Se = Det à la définition 
0j 


comme on le démontre aisément grâce à l’équation de Lagrange associée L 
ï k ” 
de l’impulsion p; associée. 


De même, c’est moins immédiat, mais on montre qu’une loi de conservation est associée à 
êL 








l'indépendance du Lagrangien par rapport au temps. En effet, supposons que 7 = 0. On peut 
alors construire la dérivée totale à partir de la différentielle, 
Rs GR COR NOR dp . dv (91) 
dt @  dtex dov dd & ? 
ce qui conduit à l’identité 
d dH 
; Le re" () 92 
(pv) = (92 


qui montre que la quantité conservée est l’hamiltonien, qui correspond dans ce cas à l'énergie totale. 
Revisitons le problème du condensateur. Dans la jauge standard, au champ électrique E est 
associé le potentiel scalaire b(x) tel que E = -®, soit 


p(x) = -Ex (93) 


où la constante d’intégration est choisie nulle. Le lagrangien vaut 





1 
L(x, à) = LE +qEx (94) 
et l'équation du mouvement 
êL d ©L 
El QE, FE Pi mi, soit mi=qE (95) 


qui est bien l’équation fondamentale de la dynamique. Du point de vue loi de conservation, on 
constate que le lagrangien ne dépend pas du temps, donc 


"08 - 
PRE 








1 
m&, soit H=pi-L ME qEx = const = 0 (96) 


avec la valeur nulle correspondant aux conditions initiales. On en déduit 


_ 2qEd 
om 


&?(d) (97) 
et c’est finalement la méthode de conservation de l’énergie mécanique que l’on retrouve ici. Notons 
que Px = _ = mi définit l'impulsion qui coïncide ici avec la quantité de mouvement et que ce 
n’est pas une quantité conservée. 

Une solution alternative consiste à faire porter le champ électrique par le potentiel vecteur A(t) 





tel que E = -dà, soit 
Alt) = -Etu, (98) 
où là encore on fixe la constante d’intégration à zéro. Le lagrangien vaut cette fois 
1 
L(x&,t) = ELLE — qEtx. (99) 
L’équation du mouvement 
êL dôL  d 
= 0 = l Et) =0, soit ï=qE 100 
ee aa ler PM à 00 
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est évidemment de nouveau l’équation fondamentale de la dynamique. Le fait que le lagrangien 
soit indépendant de x entraîne en revanche que p, est conservée, soit 





Pa = Mi + qA = mi — qEt = const = 0, (101) 


encore une fois, la valeur nulle résulte du choix de conditions. L’impulsion conservée p, diffère de 
la quantité de mouvement mt. On en d éduit 

; Et 

sd) = 7%, (102) 


m 





Notons que l’hamiltonien H = p,% — L avec à = (p; +qEt)/m et pa = 0 vaut cette fois H = =mi? 


est n’est pas une quantité conservée. C’est simplement l’énergie cinétique. 

Les diverses approches sont bien sûr compatibles entre elles, mais elles mettent l’accent sur des 
quantités conservées différentes qui sont déterminées par le choix de jauge. 
5.3.3. L’équation des ondes pour les potentiels et la jauge de Lorenz 


Pour établir des équations de propagation pour les potentiels, on part des équations avec sources 
(on pourrait les prendre nulles comme ci-dessus), 


V-E=p/e, soit V°6+@V:A = -p/e (103) 
donne l’équation de propagation du potentiel scalaire 


1 026 
V'o- De —P/Eo (104) 


à condition que soit imposée la condition de jauge de Lorenz, 
1 
V'-A+ PL 0. (105) 
De même pour le potentiel vecteur, 
ue L 2 1 1 ; 
VxB= poj + aûE, soit V'A—V(V.A)-— a V(&t6) — m0 A = —joj, (106) 


soit, en imposant la même condition de jauge de Lorenz, 


1 @2A 

VA — EE = — ho. 107 

e de H0J ( ) 

Il faut noter que les équations (103) et (106) sont bien entendu parfaitement valables, et même 

plus générales que (104) et (107) puisqu'elles ne supposent aucun choix de jauge particulier, elles 

doivent donc conduire aux mêmes résultats physiques, mais elles sont clairement plus délicates à 
manipuler et en pratique on ne le fait jamais. 


5.4. Sur l'impossibilité d’un champ uniforme dépendant du temps 


Dans l’application aux courants de Foucault, nous sommes partis d’un champ magnétique 
variable dans le temps, mais uniforme à chaque instant, 


B(t) = Bocoswtu, = B,(t)u.. (108) 


Par induction, ce champ magnétique produit une composante non uniforme du champ électrique, 
également variable, 


1 
E(r,t) = 3 Pour sinwt us = Eo(r,t) uo. (109) 
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Nous avions indiqué que ces résultats sont liés à une approximation, celle des variations temporelles 
suffisamment lentes. Pour améliorer ce modèle, nous allons considérer par exemple un champ 
magnétique créé dans une bobine parfaite parcourue par un courant variable 15 coswt, de sorte 
que Bo = uonlo si n est le nombre de spires par unité de longueur dans la bobine. Le champ 
magnétique (108) constitue l’approximation à l’ordre zéro du champ régnant dans la bobine, 


B(0 = BÜ)(r,t)u, = Bocoswtu, (110) 
et le champ électrique (109) est l’approximation d’ordre un, 
1 
E0) = EN) (r,tjuo = 3 Bowr sin wt us. (111) 
Celui-ci a été obtenu grâce à la loi de l'induction, par intégration, 
_ (rer, !)) 2 (80 (r, !)) (112) 


Cette composante d’ordre un du champ électrique donne naissance à une correction au champ 
magnétique à l’ordre deux que l’on peut estimer en utilisant la forme généralisée du théorème 
d'Ampère qui se limite ici à 
1 
VxB=- EE: (113) 


On considère une surface élémentaire rectangulaire orientée, parallèle à Oz, de hauteur h et com- 
prise entre les distances r et r + dr à l’axe Oz, dS = hdrus. Le théorème de Stokes £B - d£ — 
20, $E:4S conduit alors à 


0B% (rt) 1 0E0(r,t) 














114 
Ôr ce? ot ne 
qui donne part intégration 
(2) l w?r?2 
B£°(r,t) = — 3 Bo 2 cos wt. (115) 
À cet ordre, le champ magnétique total est alors 
1 2mr2 
B(r,t) = Bo(1-— "7 |cosutu, (116) 
4 c2 


et il n’est effectivement pas uniforme à chaque instant. On peut très bien poursuivre le calcul 
perturbatif pour produire une approximation qui sera de meilleur en meilleur aux hautes fréquences 
à mesure qu’elle comprend davantage de termes. Si l’on va à l’ordre trois pour le champ électrique, 
on applique de nouveau la loi de Faraday, 





Ô (3) _. 0 (2) 
êr ( E, (F, !)) = nr (# (, !)) ; (HIT) 
soit — 
(3) : 1 W°T® . 
E; (r, t) = ET sin Wt (118) 
et après sommation, 
1 1w?2r2\ . 
E(r,t) = 3 Bouwr (: 3 & ) sin wt Up. (119) 


De nouveau, le théorème d'Ampère généralisé fournira la correction suivante (d’odre quatre) au 
champ magnétique par résolution de 


8B% (rt) 1 0EP(r,t) 
Ôr ce? ot 
et on peut itérer assez facilement. Notons que cela engendre les termes successifs des développe- 
ments des fonctions de Bessel Ji(wr/c) et J1(wr/c) pour les amplitudes des champ magnétique et 
électrique respectivement, et que l’on obtiendrait ces solutions en résolvant directement l’équation 
des ondes. 








(120) 





16 


L2 Licence Physique Chimie 





5.5. Exercices 


Exercice 10 : Emission isotrope de charges par une sphère radioactive. 

Une sphère radioactive émet de façon isotrope des particules « (noyaux d’hélium 4He de charge 
+2/|e]) dans un milieu infini de conductivité uniforme 7. 

1.- Déterminer en tout point de l’espace la densité de courant j et le champ électrique E en 
fonction de Q(r,t), la charge électrique totale incluse dans la sphère de rayon r. 

2.- En déduire la dépendance temporelle de Q(r,t) (on notera + la conductivité du milieu), 
et exprimer j(r,t), E(r,t) et p(r,t). 

3. - En déduire la résistance électrique du volume du milieu compris entre deux sphères 
concentriques de rayons a et b. 


Exercice 11 : Cavité résonnante cylindrique. 

Un condensateur est constitué de deux armatures circulaires perpendiculaires à un axe Oz, de 
même rayon R et distantes de d. Dans tout le problème on négligera les effets de bords (R > d) 
et on supposera qu'il y a invariance par translation suivant l'axe Oz de sorte que la coordonnée z 
n'intervient explicitement dans aucune expression. 

1. - Lorsqu'une tension continue U est appliquée aux bornes du condensateur, les armatures 
prennent les charges +Q et —Q supposées uniformément réparties en surface. 

1.1 Calculer le champ électrostatique entre les armatures. 
1.2 En déduire la capacité du condensateur. 
1.3 Calculer l’énergie électrostatique emmagasinée par le condensateur. 

2.- La charge Q varie maintenant au cours du temps sous l’influence d’une tension sinusoïdale, 
ce qui donne naissance à un champ électrique variable E. 

2.1 Ecrire les équations de Maxwell et justifier l’existence d’un champ magnétique B entre 
les armatures. 

2.2 Déterminer l'orientation des champs E et B ainsi que les coordonnées dont ils dépendent. 

2.3 À basse fréquence le champ électrique peut être supposé en première approximation 
uniforme à chaque instant : E1 = Eve“*t. Calculer le champ magnétique B: engendré par la 
variation au cours du temps du flux de E; à travers une surface S: à définir. 

2.4 Calculer la correction E2 au champ électrique due à la variation temporelle du flux de 
B: à travers une seconde surface $2 à déterminer. 

2.5 Réitérer les calculs précédents pour obtenir les corrections B2 et Ez3 et en déduire un 
développement en série pour Eros. 

2.6 Tracer l’allure des approximations successives de E et B. 

3.- On donne la partie radiale de l'opérateur laplacien en coordonnées cylindriques : AE(r) = 


12 (rE(r)), ainsi que la solution f(£) sous forme de développement en série de l’équation différen- 
tielle : 
®f 1o0f 
Se = fÿ 
ae ge * /O =0, 








LE = D = 1 7 (5) + (s) _. (E) +. 


représentée ci-dessous : 

3.1 Etablir l’équation des ondes pour le champ électrique. 

3.2 Déterminer l'expression de la solution Er, t). 

3.3 On réalise une cavité résonnante en refermant le condensateur par la surface latérale d’un 
cylindre parfaitement conducteur. Déterminer en fonction de R la fréquence des deux premiers 
modes de résonance (i.e. qui correspondent à une amplitude maximale du champ électromagnétique 
dans la cavité pour des positions de la paroi cylindrique ne perturbant pas le champ électrique). 
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6. Ondes électromagnétiques 


6.1. Les équations de Maxwell dans le vide 
6.1.1. Introduction 


Les équations de Maxwell prennent une forme simplifiée dans le vide : 


(B 
rot E = E . Loi de l’induction de Faraday (121) 
divB=0, Conservation du flux magnétique (122) 
divE=0, Théorème de Gauss (123) 
1 0E 
rot B = ETS Théorème d'Ampère généralisé (124) 
c 


et c’est un fait exceptionnel que même en l’absence de sources, elles admettent des solutions prop- 
agatives E(x, t) et B(x,t) qui transportent de l’énergie, de la quantité de mouvement et du moment 
cinétique. Elles transportent la chaleur émise par le soleil, révèlent l’intérieur du corps humain ou 
facilitent les communications. Ces solutions propagatives sont souvent appelées champs libres car 
dans ce cas les lignes de champ électrique ne prennent pas naissance ni ne se referment sur des 
charges et les lignes de champ magnétique n’entourent pas des courants. Ce sont les variations 
spatio-temporelles de l’un de ces champs qui engendrent celles de l’autre et réciproquement. Ce 
résultat n’a pas été compris instantanément, car les ondes habituelles pour se propager nécessi- 
taient la présence d’un milieu dont les déformations permettaient la propagation (ondes sonores 
ou acoustiques dans un gaz ou dans un milieu condensé, mais pas dans le vide, vagues à la surface 
de l’eau, …). Rien de semblable pour les ondes électromagnétiques qui peuvent se propager dans le 
vide, c’est-à-dire sans la médiation d’un milieu matériel. 

Les équations de Maxwell sont des équations différentielles couplées pour E(x,t) et B(x,t). En 
prenant la divergence de (121) et (124) on obtient 


G(V-E)=0 et &(V-B)=0, (125) 


ces deux divergences gardent donc une valeur cosntante au cours du temps et (123) et (122) jouent 
le rôle de conditions initiales, de sorte que les champs restent à divergence nulle. 


6.1.2. Conséquences des choix de jauge 


On a vu que les champs E(x, t) et B(x,t) obéissent aux équations 


2 1®E 
Ve 0 (126) 
et ; 


Cela signifie que chaque composante cartésienne obéit à l’équation des ondes scalaire. De même 
pour les potentiels en jauge de Lorenz dans le vide, 





1 6 
2 L 
Vi  —0 (128) 
et À 
1 OA 

2 L 
Supposons, s étant un vecteur constant, que 

AL(x,t) = u(x,t}s (130) 
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(il s’agit d’une onde dite polarisée rectilignement, on discutera cette notion plus tard). L’équation 


des ondes devient D 
2 1 üU 


Cela entraîne qu’à chaque solution u(x,t) pour le potentiel vecteur de Lorenz il correspond un 
potentiel scalaire de Lorenz 


t 
bL(x, t) = —-c?s | dt'Vu(x,t). (132) 


Dans cette jauge les champs satisfont à 


E=-Vér-@Ar, B=Vx Ar. (133) 


En jauge de Coulomb on a en revanche 
V:-Ac=0. (134) 
Les équations de mouvement pour les potentiels deviennent 
V'oc = 0, (135) 
2 
1 
2 
V'Ac— de. = a V(@c). (136) 


Dans l’espace infini, la solution de l'équation de Laplace est au plus une constante et on peut la 
fixer à zéro, 


pc(x;t) = 0 (137) 
ce qui redonne l’équation de propagation ordinaire pour le potentiel vecteur de Coulomb 
2 1 Ac 
Dans cette jauge les champs satisfont à 
E=-0,Ac, B=Vx Ac. (139) 


Il est facile d’établir le théorème de Whittacker dans cette jauge : deux solutions de l’équation 
d’onde scalaire suffisent à déterminer le champ électromagnétique. En effet, le potentiel vecteur 
a trois composantes cartésiennes qui sont chacune solution de l’équation d’ondes. Or la jauge de 
Coulomb fixe une relation entre ces trois composantes et deux solutions indépendantes résolvent 
donc le problème. 


6.1.3. Ondes planes vectorielles 
Considérons l’équation générique 
1 2w 
V?w — = =0 140 
ace eu, 


et cherchons des solutions dépendant d’une seule variable spatiale (ce sont des solutions propaga- 
tives dans cette direction), w(z,t). On peut écrire 


C2w 1 d2w 
02? ce? ot? 


que l’on résoud par la méthode de d’Alembert en posant les variables indépendantes 





= (de + ce 1@) ( — c7l@) w = 0 (141) 


É=zrt+tc, m=2- ct. (142) 
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Il vient z = sin et ct = a soit Oez = à dect = 5 et dnz = à Ont = —; de sorte que 
Ô 10 1 0 
= 143 
20 da () 
0 10 1 © 
= ; 144 
On 202  2cût 4) 
L'équation (141) s’écrit simplement 
us (145) 
En 
La solution générale cherchée est donc de la forme 
w(z,t) = f(z+ ct) + g(z — ct). (146) 


La première solution se propage dans le sens de z décroissants (l'amplitude f(£) en z0 à l'instant 
to se retrouve au bout du temps t > to ultérieur en z = 20 — c(t — to) < 20) et la seconde solution 
se propage dans le sens des z croissants. 

Pour chercher E(z,t) et B(z,t) on procède comme suit. Tout d’abord 


dE, 
z 





0 (147) 


par l’équation (123) et 
0E, 
Ôt 
par l'équation (124), de sorte que E, est au plus une constante que l’on fixe à zéro car il n’y a pas 
de charges statiques dans le problème. E est par conséquent perpendiculaire à u, et l’on pose 


= Cu, -VxB=0 (148) 





E,(z,t) = fi(z+ ct), (149) 
E_(z,t) = g1(z— ct). (150) 


Le champ B(z,t) se déduit de (121), 


a on - Le = ) ru ee ss (151) 


B-_(z,i) Uuy0,91x = Ur0:91y gi(z— ct) 








Or on à ©, = de + 0h et cl © = Ôe — ©, de sorte que l’on peut relier les dérivées temporelles aux 
dérivées par rapport à z, 


a. (HO) = (AR (am ® ) 2 


0n81(n) —c7"dgi(z — ci) 





et par conséquent 


a (29) (Er) as 


Par intégration on obtient, à une fonction de z uniquement près, choisie nulle (grâce aussi à (122)) 


cB}(z,t) = -u, x fi(z+ct), k=-|klu, = kû (154) 
cB-(2,t)=u,; xgi(z-—ct) k=/klu, = kü. (155) 


Compte-tenu du sens de propagation de chacune des deux solutions, comme on l’a indiqué précédem- 
ment, on définit dans les formules ci-dessus le vecteur d'onde k pour chaque solution (nous dis- 
cuterons plus tard le module & = |k] de ce vecteur), de sorte que l’on peut finalement écrire 
l'argument de la première solution comme f, (2 + ct) = E,(—k(2+ ct)) = E,(—ku, -r — kct) = 
E.(k-r—kct) et pour la seconde g1(2—ct) = E1(k(z—ct)) = Ei(ku,-r—kct) = Ei(k-r —kct), 
soit la même forme fonctionnelle. 
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On retiendra en définitive 


k = kû (156) 
E(r,t) = Ei(k:r — kct) (157) 
cB(r,t) = à x Elr,t) (158) 


Ces ondes propagatives (E,,B.) et (E-,B-) sont appelées ondes transverses électromagnétiques 
(TEM) car elles satisfont 


E,-B,=E_.E_=0, k-E;=k-.B;=0 (159) 


et (E,B, k) forme un trièdre direct avec |[E] = c|B|. 
On appelle phase de l’onde la quantité 


p(r,t)=k:r-ckt=k:(r - cût). (160) 





La phase est constante, @(r,t) = k-ro à un instant t donné en des points r(t) tels que r(t) = ro+cût 
avec r(t) - à = ro : à + ct qui forment des plans perpendiculaires à à (pour que leur projection 
suivant à soit bien constante à t fixé) se déplaçant à la vitesse 


d 
Vo = mn = cà, (161) 


d’où la dénomination d’ondes planes. Le fait que la vitesse de phase soit égale à c = (#00) /?, la 


vitesse de propagation de la lumière, suggère bien que cette dernière est une onde électromagné- 
tique. 


6.1.4. Les propriétés mécaniques associées aux ondes électromagnétiques 
La densité d'énergie d’un champ électromagnétique quelconque s’écrit 
=" E 24 21B 2 
UM. (tt) = 5eo(lE(x 4) + c1B(x,t)°). (162) 


La densité de courant d'énergie (vecteur de Poynting S), qui satisfait à l'équation bilan 


OUE.M. 
ot 





+V.-S=-j.E (163) 


vaut 


S(x,t) = oc (E(x, t) x B(x,t)) (164) 


et elle est reliée à la densité de quantité de mouvement r(x,t) = S/c? qui, elle, obéit au bilan 


Tv. T) = —(pE + j x B) (165) 





où T est la densité de courant de quantité de mouvement qu’on appelle aussi en théorie des champs 
tenseur de contraintes. 
Dans le cas d’une onde dans le vide cela donne 





1 
ue.m.(r,t) = 3 o(Ei (kr — kct)[ + jû x Ei(k-r — kct)[?), 
= €&lEi(k:r-— kct)}, (166) 
Sr, t) LL eoclE1(k Tr — kct)|?à = UE.M.(r, t)cà, (167) 
m(x,t) = eoc |Ei(k:r — kct)?à = ugm.(r,t)c ‘à (168) 
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6.1.5. Ondes électromagnétiques planes monochromatiques 


Les ondes planes monochromatiques sont des ondes sinusoïdales, c’est-à-dire un cas particulier 
des équations (156) à (158) pour lequel la forme fonctionnelle est sinusoïdale, 


E,(r,t) = Et cos(k: r — wt + 60), (169) 
B,(r,t) = Bf cos(k-r—wt+4%0), cB° = (k/k) x Ef. (170) 


où la phase vaut w(r,t) = k:r — wt + 0 avec la pulsation 
w = ck. (171) 


On définit également une période spatiale appelée longueur d’onde À = 27/k et une période 
temporelle T = 2r/w. Dans la littérature, on introduit fréquemment la notation complexe 


Ent = cirt), £0 = Eco (rt) = R[E(r, t)] = LE +E*), (172) 
Blr,t) = B°c#@0). Bo = BO cr, B(r,t) = R[B(r,t)] = LB + B*),. (173) 


On pourrait aussi bien travailler avec les parties imaginaires qui sont également des quantité réelles. 
L'intérêt essentiel de la notation complexe est que les équations de Maxwell prennent une forme 
algébrique remarquablement simple dans ce cas, puisque V : E(r,t) = ik -Elr,t), V x E(r,t) = 
ik x Elr,t) et AE(r,t) = —iwE(r,t) de sorte qu’il vient en l’absence de charges 


) 
) 
ik x E(r,t) = iwB(r;t), ) 
ik x B(r,t) = —i(w/c)E(r, t). ) 


L'énergie transportée par l’onde plane monochromatique s'écrit en fonction des quantités com- 
plexes 





1 
UEM. — SoUREN + RIBI?) 


= Jatét (ER) + 80 -(BD)" + RE -EL + PB Beer D) (178) 





La moyenne temporelle du dernier terme est nulle de sorte qu’il reste le résultat simple 


1 
(UE.M.) = SEL (179) 


ou les crochets dénotent la moyenne temporelle, 


1 T 
co=x] de (180) 
De la même manière on exprime les moyennes 
1 x _ à 
{n) = 5e lEiPa = cKue.M>û, (181) 
{S) = cr) = Que.M.)ü. (182) 


Cela permet de définir la vitesse de propagation de l’énergie transportée par l’onde, 
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Elle coïncide avec la vitesse de phase pour l’onde plane dans le vide mais ce n’est pas toujours vrai. 
L’intensité d’une onde en un point peut être définie par 


I(x) = ÈS pat = |(S(r, 6). (184) 


Dans le cas de l’onde plane monochromatique, cela donne directement 


1 
I(r) = JE il. (185) 
C’est ce type d'expression qui fournit la raison pour laquelle en optique physique lorsque l’on 
étudie des phénomènes d’interférences ou de diffraction, on forme toujours le module au carré de 
l’amplitude complexe pour déterminer l’intensité lumineuse sur un écran ou dans un capteur. 


6.1.6. Notion de polarisation d’une onde électromagnétique plane monochromatique 


La polarisation d’une onde e.m. caractérise la direction des champs E et B dans l’espace au 
cours du temps. Par exemple pour une polarisation linéaire, ces champs conservent une direction 
constante au cours de la propagation. La polarisation d’une onde e.m. est mise à profit par exemple 
dans les dispositifs d'affichage à cristaux liquides. L’étude de la polarisation du CMB fournit des 
renseignements sur l'Univers primordial. 

On définit dans le cas général l’ellipse de polarisation. Soit les vecteurs unitaires ê1 et &2 tels 
que (ê1, &2, à) soit direct. Ils forment une base pour Ef, soit pour l’onde monochromatique, 


E(r,t) = (E1ê81 HÉnes entrer) (186) 
et on peut poser | | 
E1= Ajet, Eo = Ape®, p=k.r-uwt, (187) 
de sorte que 
R[E] = À; cos(ç + 01)81 + A) cos( + 02)@, 
= E: (r, t)è: + Es (r, t)82. (188) 


Avec un peu de trigonométrie, on montre que 


cos(p + 01) sin 02 — cos(çp + 02) sin 01 = coswsin(02 — 01), (189) 





cos(( + 1) cos d2 — cos(ç + 02) cos 01 = sin y cos(02 — 01), (190) 


soit en sommant les carrés de ces deux égalités 


FE: : FE) 2 E: E2 :. 2 
(à) L (à) er À) SU (4 





qui est l’équation d’une ellipse dans le plan Æ1,E2 dont le grand axe est incliné par rapport à E: 
d’un angle a tel que 





2A1A2 
tan 2a = cos(02 — 01). (192) 
A? — A 
Dans le cas où 62—01 = mn, m=0,+1...,les composantes orthogonales du champ électrique 


sont en phase ou en opposition de phase, l’ellipse dégénère en ligne droite et cela correspond à une 
polarisation linéaire car Er, t) pointe dans une direction fixe Vr,t et l’on peut écrire 


R[EE] = (Arê1 + Avêo) cos(k : r — wt + 6). (193) 


Si en revanche 02 — 61 = mr/2, m = +1,+3..., les composantes orthogonales du champ 
électrique ont même amplitude mais sont en quadrature de phase et on obtient une polarisation 
circulaire. 

Les ondes polarisées rectilignement, tout comme les ondes polarisées circulairement constituent 
des bases pour la représentation générale d’une onde. 
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6.1.7. Ondes cylindriques dans le vide 


Nous avons traité uniquement jusqu’à présent le cas des ondes planes. Ce sont des solutions 
particulières des équations de Maxwell dans le vide, mais il en existe de nombreuses autres ayant 
des symétries différentes imposées par exemple par des conditions aux limites déterminées par les 
conditions de production de l’onde. 

On considère ici des solutions de symétrie cylindrique. La résolution exacte dans ce cas est 
difficile. Elle fait appel à des fonctions spéciales de la physique mathématique (les fonctions de 
Bessel) que nous n’aborderons pas ici. Nous allons en revanche traiter le problème de manière 
approchée. 

On cherchera donc une solution réelle des équations de Maxwell de la forme 


E(r,t) = Eof(r) cos(kr —wt)u,, f(r) ER (194) 


en coordonnées (r,0,z). C’est une onde monochromatique polarisée rectilignement. On se place 
dans le vide si bien que w = ck. L’onde se propage de manière radiale depuis l’axe Oz vers 
l'extérieur, elle est invariante par rotation d’angle 0 autour de cet axe et par translation le long de 
Oz, les surface équiphase (surfaces d'onde) sont des cylindres d’axe Oz. La fonction f(r) est sans 
dimension. 

Pour exprimer B(r,t) on utilise l'équation de Maxwell-Faraday 




















VXxE = -0,B (195) 
avec pour le rotationnel 
10E, 0E, 
VxE 20 Ur 7 (196) 
où seul le deuxième terme subsiste et d’où l’on déduit 
0B; 
— = 0, (197) 
ot 
°0B, 0E, d 
8 sp UX cos(kr — wt) — kf(r) sin(kr — wt) |, (198) 
ot Ôr dr 
0B, 
= (0. 1 
ot (199) 
Par intégration sur t et en l’absence de champ magnétique constant on obtient 
B(r,t) = Bo(r,t)u, (200) 
Eo [à 
Bo(r,t) : ë sin(kr — wt) + kf(r) cos(kr — ur) (201) 
w r 


Pour fermer le système d'équations (et obtenir une contrainte sur la fonction f(r)), on utilise 
l'équation de Maxwell-Ampère, V x B = LOE. 








10 10 
A 02) 
r Ôr r Oz 
ce qui conduit à l’équation 
1 0 1 0 1 °0Bo 1 0 
pie E, soit -B = E.. 203 
r Ôr te) ce? Ôt re | Ôr c2 Ot Ge 


En remplaçant les champs par leurs expressions en fonction de la fonction f(r) on parvient à 
l'équation suivante 


1[df . dia df 
. Ë sin(kr — wt) + T2 sin(kr — wt) + kr cos(kr — wt) 
df 2 ; Le . 
+kf(r) cos(kr — wt) + kr cos — k°rf(r)sin(kr —wt)| = — sin(kr — wt). (204) 
r c 
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La relation de dispersion permet de simplifier le second membre avec le dernier terme du premier 
membre et si l’on sépare les coefficients de cos de ceux de sin on obtient deux équations, 


1 

FPE POS O (205) 

2kr f'(r) + kf(r) = 0, (206) 
où il est préférable de tout diviser par k pour rendre explicite les simplifications dans la limite de 
grande distance kr » 1, soit 

ES") + Pr) = 0 (207) 

k kr . 

2r f'(r) + f(r) = 0. (208) 


Dans cette limite kr >» 1, le deuxième terme de la première équation est négligeable, ce qui impose 
que le premier le soit aussi. On peut alors chercher une solution en loi de puissance de la seconde 
équation, f(r) + ar* qui conduit à la solution f(r) — 1/4/r de sorte que dans cette approximation 
il vient les solutions approchées pour l’onde cylindrique 





1 
E(r,t) = NE. cos(kr — wt)u,, (209) 
1 


VT c 


On remarque qu'avec f(r) + 1/4/r, les deux termes de la première équation sont bien négligeables 
à grande distance. 

Le même résultat approché peut être obtenu au moyen de considérations énergétiques. A l’aide 
des expressions (194) et (201), on peut calculer le vecteur de Poynting 


B(r,t) 








cos(kr — wt)u. (210) 


LE(r, t) x B(r,t) 
Ho 


1 
= ——E,(r,t)Bo(r,t) 
H0 


S(r,t) 


L 1 d df 
= Le 20) cos? (kr — wt) mi = 





sin(kr — wt) cos(kr — «| u,. (211) 
En prenant les moyennes temporelles (cos?) = + et (sin cos) = 0 il vient 


S0) = Pr) (19) 


On peut alors calculer la puissance émise à travers un cylindre de hauteur unité et de rayon r, 


1 27 
Puissance = | de | rd0<S) - u, 
0 0 


Eé 
mec (r) 


= eocrnESrf”(r). (213) 
En l’absence de phénomènes dissipatifs, cette puissance est égale à une constante P5, d’où 


f(r)= _ (214) 


TEocEër 


qui confirme l’expression précédente. 
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6.1.8. Exercices 


Exercice 12 : Station radiophonique. 

Une station radiophonique émet sur 1829 m (grandes ondes). Supposons qu’elle envoie une onde 
dans la direction Ox, polarisée suivant Oy, de champ électrique Eg = 3.10% V/m dans la région 
où se trouve le récepteur. 

1.- Calculer w et B. 

2.- Comment doit-on placer un récepteur constitué d’un cadre rectangulaire de surface S = 
100 cm? comprenant N — 10 spires, pour recevoir la fem e maximale ? Calculer e dans ces 
conditions. 

3.- Calculer le vecteur de Poynting et la puissance émettrice sachant que le récepteur est situé 
à la distance À = 100 km de l’émetteur et que la station émet sensiblement la même énergie dans 
toutes les directions. 


Exercice 13 : Onde électromagnétique captée par une antenne. 

Une onde électromagnétique plane polarisée rectilignement est captée par un cadre rectangulaire 
de côté a comportant n tours de fil. 

1. - Comment orienter le cadre et choisir a pour que la fem e soit maximale ? 

2.- Comment déduire de la valeur efficace de la fem l’amplitude du champ électrique ? On 
note Æ cette valeur. L’émetteur est à une distance D »5 a du cadre. Calculer sa puissance P. 

AN :v = 100 MHz,n = 10, E = 0.2 mV, D = 50 km. 


Exercice 14 : Equations algébriques. 
Soit un champ de vecteurs 
E(r, t) = Epe(FT-wt) 


où le vecteur d'onde k à pour composantes k,, k,, k. et le vecteur (constant avec F et t) E, Ex, 
Eoy; Eo-. Exprimer div E et rot E en fonction de k et E. 


Exercice 15 : Dipôle sinusoïdal. 
Un dipôle sinusoïdal oscille dans la direction Oz autour du point O : 


U(t) = acoswt = Re (ae!) 


1. - Calculer en coordonnées polaires les composantes du potentiel vecteur et le potentiel 
scalaire. 
2.- En déduire les champs électrique et magnétique et le vecteur de Poynting. 


6.2. Quelques notions sur la propagation dans un milieu matériel ou en son voisinage 
6.2.1. Relations constitutives dans les milieux 


L'étude de la propagation dans les milieux matériels dépasse largement le cadre de ce cours, 
mais quelques résultats peuvent être obtenus par l'intermédiaire d'exercices guidés. Les milieux 
les plus simples sont linéaires, homogènes et isotropes. On distingue les milieux transparents ou 
diélectriques, les conducteurs (métaux ou par exemple supraconducteurs), les milieux magnétiques 
simples (paramagnétiques). Bien entendu de nombreux matériaux n’entrent pas dans ce cadre, par 
exemple les ferromagnétiques. 

La présence d’un milieu peut aussi affecter la propagation de l’onde dans le vide en son voisinage. 
Cela conduit à des résultats intéressants que l’on met à profit dans les guides d’onde par exemple. 

On notera que dans les cas les plus simples, un milieu transparent sera caractérisé par une 
permittivité diélectrique 

€ = ErEo 


distincte de celle du vide et un milieu magnétique par une perméabilité magnétique 


H— HrHo 
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également distincte de celle du vide. 
Dans un conducteur, une relation importante est la loi d’'Ohm locale 


j=0E 


et dans d’autres matériaux, comme les supraconducteurs, cette relation doit être modifiée pour 
rendre compte de la présence de courants non dissipatifs, 


j = -aA. 


6.2.2. Effet de peau 


Dans un métal, en raison de la grande conductivité, le champ électrique ne pénètre que très 
peu. On appelle ce phénomène l’effet de peau. Il dépend de la fréquence comme on va le voir, 
et plus la fréquence de l’onde électromagnétique est élevée, plus la profondeur de pénétration du 
champ électrique est petite. 

On considère pour cela un métal de conductivité o occupant le demi-espace x > 0 et pour 
lequel on cherche une solution des équations de Maxwell sous la forme d’ondes planes progressives 
polarisées de pulsation w. Dans un métal, dans certaines gammes de fréquences, le courant de 
déplacement €o est négligeable devant le courant de conduction j = oE. On supposera que c’est 
le cas ici et l’on quantifiera les gammes de fréquence en question un peu plus tard. 

On cherche une expression de la forme : 


E(x,t) = Eof(x) cos(kx — wt)u, 
ou, en notation complexe : 
E(at) = Eof(x)e 7 u, 


où u, est le vecteur unitaire sur Oz et f(x) est une fonction réelle que l’on cherche à préciser. 
L’onde se propage suivant Ox, elle est polarisée sujivant Oz, on peut donc chercher au moins une 
composante de B suivant Oy. 

En coordonnées catrésiennes, l’équation de Maxwell-Faraday conduit à 


0E,(x,t)  OB,(x,t) 





Ôx ot 
que l’on intègre en 
1 i(kx—wt) df : 
B(x,t) = R|-—Ee — +ikf]|u, 
iw dx 
L’équation de Maxwell-Ampère, compte-tenu des hypothèses indiquées plus haut, prend la forme 
VxB= /H100E. 
Elle conduit ici à 0B, («0 
Lit 
— = L00E, (xt), 


soit, avec la forme précédente de B,(x,t) et de E,(x,t) l'équation diffférentielle pour f(x) 








Loow f (x) 2 +i (TS 7) , 


dx? 


On identifie les parties réelle et imaginaire pour obtenir deux équations, 








por f(e) = 28 7), 
0 = 0 — K2 f(x). (215) 
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La seconde équation admet des solutions sous la forme d’exponentielles réelles dont on exclut la 
solution divergente et il vient 


f(x) = AeTFe 
que l’on reporte dans la première pour identifier 


PE, 5=-x! 

2 ” ‘ 
La constante Ô représente l’ordre de grandeur de la distance sur laquelle l’onde électromagnétique 
peut pénétrer dans le métal. Aux fréquences v = 102, 10%, 104 et 10°Hz elle vaut respectivement 


6.6, 2.1, 0.66 et 0.21 mm. 


6.2.3. Exercices 


Exercice 16 : Effet de peau dans un conducteur. 

Les champs E et B à l’intérieur d’un conducteur de conductivité 7 varient dans le temps. 

1. - Quelle est l'équation différentielle vérifiée a le vecteur de densité de courant j quand 
on néglige la densité de courant de PR . Ep a - Pour quelles fréquences, l’équation obtenue 
est-elle valable ? (On donne pour le cuivre + = 6-: “107 (QG :m) 1.) 

: 2. - En utilisant la notation complexe des grandeurs électromagnétiques, on pose : j = 
3() exp(—iwt) en supposant que l’amplitude complexe du vecteur densité de courant ne dépend 
que de la coordonnée z. 


Quelle est la forme de j(z) ? On introduira la profondeur de peau : Ô = 





ne Indiquer la 
signification physique des différents termes. 
3.- Evaluer 6 dans le cas du cuivre pour des fréquences de 50 H2, de 1 MHz et de 10* MHz. 


Conclusions ? 


Exercice 17 : Supraconducteur : équation de London et effet Meissner. 
Un matériau supraconducteur présente les constantes électrique et magnétique du vide (£0 et 
Lo), mais la loi d'Ohm y est remplacée par la relation constitutive dite de London : 


G est la densité de courant, À le potentiel vecteur, À une constante caractéristique du matériau 
étudié). 

1.- Trouver l’équation différentielle satisfaite par le champ magnétique B. 

2.- Quand B ne dépend pas du temps, comment variera-t-il au voisinage de la surface (supposée 
plane et infinie) d’un supraconducteur immergé dans un champ extérieur parallèle à cette surface 
? 

3.- En résine variable, la densité de courant j est somme de deux termes : l’un, a ant 
à loi d'Ohm (j 11 = oË où o sera supposée indépendante de la fréquence) et l’autre, FLE obéissant 
à l’équation de London. Donner la relation de dispersion pour une onde plane de DHÉaG on w et 
nombre d’onde k propagée dans le milieu. 


Exercice 18 : Propagation d’une onde électromagnétique parallèlement à un plan conducteur, 
pression de radiation. 

Le demi-espace z < 0 est occupé par un conducteur parfait, le reste étant occupé par le vide. 
On étudie la progation d’une onde électromagnétique plane dans la direction Ox. 

1. - Exprimer les champs E et B. 

2.- Calculer les densités de charges et de courants superficiels. 

3. - Quelle est la pression de radiation exercée par l’onde sur la surface métallique ? 
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Exercice 19 : Propagation d’une onde électromagnétique entre deux plans conducteurs. 
Une onde électromagnétique se propage dans le vide entre deux plans métalliques situés en 
z = 0 et z = a. L’onde étant transverse électrique, le champ électrique correspondant s’écrit : 


E = Eosin (=) cos (kx — wt)j 

1. - Etablir l’équation de propagation du champ électrique puis à l’aide de la forme [1] en 
déduire la relation de dispersion k = f(w). Tracer la courbe k = f(w). 

2.- En généralisant la relation À = c/v, valable pour la propagation dans l’espace vide infini, 
déterminer l’indice de réfraction ”virtuel” de ce milieu dispersif (4.e. pour lequel l’indice est fonction 
de la fréquence). 

3.- Calculer le champ magnétique associé à cette onde. L’onde est-elle transverse magnétique? 

4.- Déterminer le vecteur de Poynting T. 

5.- Calculer la densité d'énergie associée aux champs et en déduire la vitesse de propagation 
de l’énergie. Conclusions. 


7. Formes différentielles et équations de Maxwell 


7.1. Quelques notions mathématiques 
7.1.1. 1—, 2-—, et 3—formes dans R° 


Les formes différentielles constituent un langage mathématique particulièrement adapté à l’élec- 
tromagnétisme. Nous en présentons les bases dans cette section, même si cela n’est en aucun cas 
exigible au niveau L2. 

On se placera dans le formalisme à 3 + 1 dimensions dans lequel l’espace M est une variété à 
D = 3 dimensions supposée avoir des coordonnées {x/};_123. On n’expliquera pas pourquoi les 
coordonnées sont notées avec un indice “en haut”, c’est lié aux notations tensorielles, mais nous 
allons l’admettre et travailler de cette manière. 

On appelera forme différentielle de degré p ou plus simplement p—forme la quantité notée & 
qui, intégrée sur un domaine p—dimensionnel . (un domaine de M), délivre une quantité scalaire 


8, (Notez la position de p pour la p—-forme et pour le domaine d'intégration. Ce n’est pas une 
notation très répandue, mais elle simplifie la compréhension) 


[Le = $. (216) 


P 


Une 1-forme w integrée le long d’une ligne £ = ( est équivalente à la circulation le long d’une 


1 
Jo- [ua (217) 
Q £ 


avec w = we) un vecteur de l’espace 3-dimensionnel cotangent (là encore on n’explique pas la 
raison de cette notation mais on admet d’un co-vecteur a des composantes dénotées à l’aide d’un 
indice “en bas” et un vecteur a des composantes avec un indice “en haut”) et dx = drey l'élément 
de longueur de l’espace tangent. À trois dimensions, une 1-forme peut ainsi s’écrire 


ligne d’un champ de vecteurs 


& = widr! + wodx? + w3dx° (218) 
et elle correspond à la circulation infinitésimale du co-vecteur 
ü 2 w-dx. (219) 


Ici, la notation = signifie “correspond à” et nous la prenons au sens d’une correspondance exacte, 
ie. une égalité, comme dans (218). 
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sq : 2 : : 1 
Une forme est indépendente du choix de coordonnées. L’invariance de w sous le changement de 
42 j j! 1 
coordonnées x — x/ montre que w; se transforme comme un tenseur covariant de rang 1, 


0x 


Oxi 


W j ou W. (220) 


DRE one : 
Une 2-forme et une 3-forme, w et w, intégrées respectivement sur une surface S = Q, ou sur un 
2 


volume Y = (, sont équivalentes à des intégrales de surface, 
3 


Jé-fus (221) 
Q S 
et de volume, 


1 & = Î wdV. (222) 


Les 2- et 3-formes peuvent s’écrire comme 


2 
w = widr! dr? + wosdr? dr + ws1dx° dx", (223) 
3 

w = wissdx dr? dx*, (224) 
mais dans ces deux expressions, une orientation doit être choisie pour spécifier complètement le 
signe des intégrales. Par exemple le membre de droite de (221) suppose une orientation pour 
dS. C’est un example d’intégrale bidimensionnelle orientée qui se transforme par changement de 
variable selon 





= Hé no ED do 
12 fewdrdu = | F(a(u, v), y(u, tu y d (225) 


S(u,v) 
où le déterminant jacobien n’est pas pris en valeur absolue précisément en raison de la spécificication 
d'orientation. On peut observer par exemple avec u = y et v = x qu’il vient | drdy = —\dydx 
et cela motive l'introduction du produit extérieur À (aussi appelé produit wedge) de sorte que la 
propriété d’antisymétrie soit satisfaite, 
da? À dx = —dx* À dx?, et en particulier dx? À dx? = 0. (226) 


Une 2-forme peut ainsi être écrite proprement à l’aide de 


wiodx! dx? — Que A dr? + und? A dt}, War = —wi2 (227) 
comme — 
w = 7 ik da? A dx, (228) 
et pour une 3—forme, 
= Suundr À di À dr”. (229) 
Leurs composantes se transforment comme des tenseurs covariants de rang p, 
j 5h 
Comme pour la 1-forme, on suppose que la correspondance 
à 2 w.dS (231) 
ÿ 2 wdV (232) 


est une égalité. 
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7.1.2. Produit de Hodge et dérivée extérieure 


Plusieurs opérations entre formes différentielles peuvent être réalisées, et parmi elles, l’opérateur 
de dualité Hodge (ou produit étoilé) et la différentielle extérieure (ou dérivée extérieure) jouent un 
rôle essentiel. 

La différentielle extérieure d augmente le rang d’une forme d’une unité et pour une forme définie 
par 





1 
= —un.….pdtt A+: À da”, (233) 
DE 
on aura Ï à 
de == Pl grd A dx! À: A dx?. (234) 
pl  Oxi 


L'opérateur d agissant sur un produit obéit à une loi de Leibniz généralisée, 
p1 p 1 p À 
d(wb) = (dw)0 + (—-1)Pwd6. (235) 
Il est facile de voir que dé produit un gradient, du produit un rotationnel, et dé produit une 


divergence. 


Exercice 20 : Vérifier la proposition ci-dessus dans un système de coordonnées curviligne, comme 
les coordonnées sphériques. 


0 
Pour une 1—forme ® = ®, on a 


; 0 20 0 





Ô 

_ d (A 2 
do Gr r+ sd or (236) 
grad  - dd une 1—forme, d£ est un déplacement général (237) 


= (gradé:u,)dr + (grad à - us)rdô + (grad @ : u,,)r sin Ode 


soit par identification 





_ 0@ 1 09 1 0 
rar or Fe 00  ‘ rsinÿ op Ke C8) 
1 
Pour une 1—forme E, on a 
1 
E = E,dr + Eodô + E,de (239) 
= E,dr + Eprd0 + Er sin /dp (240) 





où les composantes “bar” sont celles du champ (électrique ici) physique, c’est-à-dire avec les facteurs 
métriques et les mêmes unités pour chacune des trois composantes, ce qui n’est pas le cas des 
composantes nues. On forme la dérivée extérieure, 
1 ÔE; QE, 
dE — T0 À dr + 
00 0p 
ÔE, 0E, 
+ dr À d0 + dep À dô 
Ôr 0 
0E, 0E, 
Ôr 00 


= rotE:dS une 2-forme, dS est une surface générale (242) 
= (rotE:u,)rdôr sin Ode + (rot E - uo)r sin Odçpdr + (rot E - u,,)drrdô (243) 





dp À dr 








dr À dp + dô À dp (241) 


et par identification par exemple 


> 
—— — E;. 
r sin 0 0 " 





1 0E E, 
(rot E : u,) = (£ Gi 0e 


INC 
30 2e ) (sMm0E,,) — 


r? sin 0 r sin 0 00 
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On procède de même pour les autres composantes. 


2 
Dans le cas d’une 2—forme B on écrit 


2 
B =  B;odr À dô + Bo,0 À dy + Bd À dr (245) 
= B,rdôr sin 0dp + Bor sin Odpdr + B,,drrd0 (246) 


et de nouveau les composantes B; sont celles du champ physique, liées à celles de la 2—forme, par 


exemple Bo, = r?sin 0B,. On exprime maintenant dB, 











2 B; B Br 
LE: TT dd 0 0 ie 87 0 À de À dr 
0p — — Ôr 0) 
dr A dO A do 
= es + (r?gin@ sn 08 )| dr À dd À d (247) 
[op 7 ôr "7 0 t F 
= divBdV une 3-forme (248) 
= divBr°sin 0drdôdp (249) 
L'expression de la divergence s’ensuit : 
1 à - 1 à = 1 0B 
divB = ——(r28B,) + ———(sin0B +. 2 
Fe r? Ôr VAPREE r sin O 00 GÉFBUN rsim0 0 C0) 


Le produit de Hodge (dénoté x3) agissant sur une p-forme & (p < 3 ici) délivre une (3—p)-forme 


3—p re ns 
w . Dans un espace métrique de tenseur métrique 


9jk — €j ° €k; gi5g Æ 6}, (251) 
on définit 
1-forme RE w;dx° (252) 
ie = à — order” A dx* 
1 , 
= 3 Vis g"ur dr? A da*, (253) 
2 1 ï £ 
2-forme WU — ji is dr A dx? 
430 = & = ydr* 
= V9Ei5r9 "9 "wimdr", (254) 
3 1 5 ; k 
3-forme wW — grise de A dx? À dx 
3 0 
X3W — YU 
= VGEi5kg" 9 9 "uwimn- (255) 


On voit apparaître des facteurs métriques dans la définition du Hodge, c’est donc une opération 
qui intervient après qu’on à introduit une structure métrique sur la variété M. C’est un point 
technique qui ne pose pas de problème ici car pour les physiciens, les mesures de longueurs (ou de 
durées) font partie des opérations de base et la structure métrique est généralement installée. 
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7.1.3. Lemme de Poincaré et théorème de Stokes 
Un résultat important est le lemme de Poincaré d(dt) = 0 pour toute forme, ou encore 


d = 0. (256) 


p 
On en déduit que si une forme est fermée, i.e. obéit à dû = 0, alors elle peut être écrite comme 


p #1 _1 
0 = d(fw ) et elle est dite exacte (on peut appeler lw un potentiel, puisqu'elle est elle-même définie 
à la dérivée extérieure d’une (p — 2)-forme près). Un cas particulier évident est celui où p = D, 


de = 0. (257) 
Enfin, le théorème de Stokes prend une forme très compacte en termes de formes différentielles, 
de = | & (258) 
Q Ce) 
p+1l p+1l 


où © {À est le bord de la sous-variété (p + 1)-dimensionnelle Q. 
p+1 p+1 


7.2. Application aux équations de Maxwell 
7.2.1. Les quatre équations fondamentales 


Le théorème de Stokes joue un rôle central dans l’élaboration de l’électromagnétisme. Nous 
allons suivre ici une démarche élaborée par Hehl et Obukhov. 

La charge totale Q comprise dans un domaine de volume V de l’espace est un scalaire et peut 
être écrite comme l'intégrale de volume d’une densité de charges p, 


Q= [nv (259) 


d’où l’on déduit que cette dernière est une 3-forme (dans l’espace ordinaire à 3 dimensions), D telle 
que 


3 
Q= | à. (260) 
v 
S'agissant d’une 3-forme dans l’espace ordinaire à 3 dimensions, il en résulte que 


dp=0 (261) 
3 2 
ce qui demande que p soit la dérivée extérieure d’une 2-forme que nous notons D, 


2 
dD = p. (262) 


2 
Ce n’est rien d’autre que l’équation de Maxwell-Gauss et la notation D suggère un champ électrique 
auxiliaire (aussi appelé déplacement électrique). 


On peut aussi considérer les variations temporelles de Q. La dérivée temporelle de est non 
nulle si des charges entrent dans le volume V ou le quittent, une propriété de conservation que l’on 
écrit d 

np rav--fius (263) 
dt Jy 
S 


où S = OV est le bord du volume considéré et j est la densité de courant. L'intégrale de surface (le 
flux) de j à travers une surface quelconque X est l'intensité du courant (orientée) 7 (un scalaire) 


2 
traversant la surface, d’où il découle qu’une 2-forme 7 peut être définie via 


Te Î 3. (264) 
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À laide du théorème de Stokes, l’équation de conservation de la charge éléctrique (263) devient 


3 2 3 2 
[ a+ | = [ (+ a —0 (265) 
v ov v 
et comme le domaine V est arbitraire, 
dp + dj = 0. (266) 
2 
En combinant cette équation et la loi de Gauss (262) on voit que 5 + ©.D est une 2-forme exacte, 
9 2 
d(3 + &D) = 0, (267) 
qui permet d'introduire une 1-forme comme un potentiel , 
1. 2 2 
dH = + 0,D. (268) 
1 
C’est la loi de Maxwell-Ampère avec la notation H pour le champ auxiliaire magnétique (ou 
excitation magnétique). 
Les deux équations de Maxwell restantes peuvent être obtenues en suivant un raisonnement 


analogue, mais avec le flux magnétique et la différence de potentiel comme quantités scalaires de 
départ. En effet, 


x 
Il 


1 B-dS, (269) 
2 


V = JE. (270) 
£ 


2 1 
montrent qu’une 2-forme B et une 1-forme E peuvent être définies de sorte que 


2 
& — Je (271) 
X 
1 
V = |lE. (272) 
£ 


On peut supposer qu’une densité linéique de courant de flux j°? existe pour la variation de flux 
magnétique à travers une surface fermée comme une loi de conservation, 


a] 8-4s-- } j° - dx, (273) 
del 
L=0S 


ou, en termes de formes différentielles, 


2 2 
| B + | 5° = | (dB + dj?) = 0. (274) 
S OS S 


La densité linéique de courant de flux est une 1-forme, comme le champ électrique et elle a les 
mêmes dimensions que celles du champ électrique. Ayant les mêmes propriétés géométriques et 
dimensionnelles, on peut supposer qu’il s’agit de la même quantité, à un facteur numérique multi- 
plicatif sans dimension près (qui se révèle valoir 1), On obtient ainsi 


2 1 
dB + dE =0, (275) 


que l’on reconnaît comme la loi de l'induction de Maxwell-Faraday. 
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On applique maintenant l’opérateur d à cette équation, de sorte qu’en utilisant le lemme de 
2 
Poincaré il vient (dB) = 0. La constante d'intégration doit s’annuler, sinon il existerait des 


1 

53 à PE : À : k : 1 « NÉ 
charges magnétiques, les monopôles magnétiques, qui empêcheraient d'identifier 9% à Æ, mais il 
n'existe aucune trace expérimentale de tels monopôles. On arrive ainsi à la dernière équation de 
Maxwell, 


2 
dB = 0. (276) 


Un point très important à mentionner est que les équations (262), (268), (275) et (276) sont très 
générales. Elles ont été obtenues sur la base d'arguments géométriques et de lois de conservation 
et sont écrites dans un langage indépendant des coordonnées, elles admettent donc des formes 
universelles valables dans n'importe quelle variété tridimensionnelle (y compris non euclidienne). 
Cette présentation est dite pré-métrique, car aucune structure métrique n’a été supposée. Pour 
être plus concret sur ce point, supposons que les trois coordonnées sur M s’appellent p, 4, z (cela 
suggère évidemment les coordonnées cylindriques, mais à ce niveau, z pourrait très bien être un 
angle !). L’équation (262), avec 

















2 
D = D,dp À dp + D,zdp À dz + D,,dz À dp (277) 
et . 
p = gdp À dy À dz, (278) 
conduit à Se mn 35 
pe ?£ RÉ = 279 
Ôz u Op L Ce) ’ (AS) 
Avec 
1 
H = H,dp + H,dç + H,dz (280) 
et é 
31= odp À dp + ipzdp À dr + i,9d2 À dp, (281) 
(268) correspond à trois équations, 
0H, 0H, . 
SP}; 282 
ôp êp Lys ( ) 
CH, CH, 
on — Lzp) 2 
êz op lzp ( 83) 
0H, OH, . 
an 284 
êp êz lo ( 8 ) 


De même, (275) fournit trois équations 


0, 0E By 

















: 2 
D en 0 (285) 
CE: 0B2  0BEr- 
Ôz Ôôp ôt % (266) 
0Ë,, 0Ë... 0Bis 
01) ôz  ôt ® (87) 
et enfin (276) s’écrit 

Be + Be: + °B2 — 0: (288) 

CA 0p Ôp 


Aucune de ces équations ne contient de facteur métrique, les trois coordonnées sont totalement 
arbitraires et pourraient être appelées e, x, D ! 
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7.2.2. Nécessité d’une structure métrique 


1 
Au niveau pré-métrique, il y a donc quatre champs fondamentaux décrits par les 1-formes H et 

1 2 2 

E, et par les 2-formes D et B, et a priori aucune relation entre elles. Le système d'équations n’est 


donc pas fermé et il manque des informations. Les relations entre les champs fondamentaux (autres 
que les équations de Maxwell) apparaissent dès lors que l’on introduit les relations constitutives 
du milieu dans lequel règnent les champs, par exemple, en supposant la linéarité, essentielle à 
l’électromagnétisme et assise sur des bases expérimentales solides, 


D = «E, (289) 
H = y !B. (290) 


Ces deux relations toutefois ne sont pas tout-à-fait correctes au sens des formes différentielles et 
doivent être traduites en 

2 1 
D = E(x3E), (291) 


1 2 
H = px), (292) 


où le produit de Hodge est nécessaire pour restaurer des formes de même degré de part et d’autre 
des égalités. Comme ces relations impliquent le produit de Hodge, elles présupposent une métrique 
et elles dépendent par conséquent du choix du système de coordonnées. De ce point de vue, même 
le vide est une structure et induit une relation entre les quatre champs (avec € = €0 et du = Ho). 


Exercice 21 : Cylindre chargé. 


On considère un cylindre de longueur infinie et de rayon a portant une densité de charge 
2 1 
uniforme. Déterminez les composantes des formes D et E dans l’espace euclidien en coordonnées 


cylindriques. 


Exercice 22 : Cylindre parcouru par un courant uniforme. 
Dans la même géométrie et pour un cylindre parcouru par un courant uniforme dans la direction 
1 2 
de son axe, déterminez les composantes des formes H et B 
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